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1. Introduction 



We consider scalar conservation laws of the form 



(1.1) 



dtu + dxf{u) = 0, u = u{x,t) eR,t>0, 



where the flux-function / : M ^ M is a given, smooth mapping. It is well- 
known that initially smooth solutions of (1.1) develop singularities in finite 
time and that weak solutions satisfying (1.1) in the sense of distributions 
together with a suitable entropy condition must be sought. For instance, 
when the initial data have bounded variation, the Cauchy problem for (1.1) 
admits a unique entropy solution in the class of bounded functions with 
bounded variation. (See, for instance, [8].) In the present paper, we are 
primarily interested in shock waves of (1.1), i.e. step-functions propagating 
at constant speed. 

Entropy solutions of (1.1) can be obtained as limits of diffusion or relax- 
ation models. For instance, under the sub-characteristic condition [9] 



and when the relaxation parameter e > tends to zero it is not difficult to 
check that solutions of 



converge toward entropy solutions of (1.1). More precisely, the first compo- 
nent u := lim^^o^^^ is an entropy solution of (1.1) and f{u) := lim^^o^^^ is 
the corresponding flux. See, for instance, Natalini [11] and the references 
therein for a review and references. 

The Chapman-Enskog approach [2] allows one to approximate (to "first- 
order") the relaxation model (1.3) by a diffusion equation ((1.4) below). 
More generally, it provides a natural connection between the kinetic de- 
scription of gas dynamics and the macroscopic description of continuum 
mechanics. The Chapman-Enskog expansion and its variants have received 
a lot of attention, from many different perspectives. For recent works on 
relaxation models like (1.3), Chapman-Enskog expansions, and related mat- 
ters we refer to Liu [9j, Caflisch and Liu pj, Szepessy [13j, Natalini [llj . 
Mascia and Natalini [lOj, Slemrod [12], Jin and Slemrod [6], Klingenberg 
and al. [7], and the many references therein. 

Our goal in this paper is to initiate the investigation of the validity of the 
Chapman-Enskog expansion for discontinuous solutions containing shock 
waves. This expansion is described in the literature for solutions which are 
sufficiently smooth, and it is not a priori clear that such a formal procedure 
could still be valid for discontinuous solutions. This issue does not seem 
to have received the attention it deserves, however. Note first that, by the 





(1.3) 



dtUe + d^Ve = 0, 

dtVe + dxUe = - 



{f{Ue) - V, 
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second equation in (1.3), we formally have 



f{ue) - e {dtv^ + dxUe) 

f{u,) - e {-f'{u,) d^f{u,) + a" d^u,) + 0{e^) 



as long as second-order derivatives of the solution remain uniformly bounded 
in e. Keeping first-order terms only, we arrive at the diffusion equation 



This expansion can be continued at higher-order to provide, for smooth 
solutions of (1.3), an approximation with higher accuracy. When solutions 
of (1.3) cease to be smooth and the gradient dxU^ becomes large, the terms 
collected in O(e^) above are clearly no longer negligible in a neighborhood 
of jumps. The validity of the first-order approximation (1.4), as well as 
higher-order expansions in powers of e, becomes questionable. 

The present paper is motivated by earlier results by Goodman and Majda 
[3] (validity of the equivalent equation associated with a difference scheme) , 
Hou and LeFloch [5] (difference schemes in nonconservative form) , and Hayes 
and LeFloch [i] (diffusive-dispersive schemes to compute nonclassical en- 
tropy solutions). In these three papers, the validity of an asymptotic method 
is investigated for discontinuous solutions, by restricting attention to shock 
waves with sufficiently small strength. This is the point of view we will 
adopt and, in the present paper, we provide a rigorous justification of the 
validity of the Chapman-Enskog expansion for solutions containing shocks 
with small strength. 

Specifically, restricting attention to traveling wave solutions of the re- 
laxation model (1.3), the first-order approximation (1.4), and the associ- 
ated second-order approximation (see Section 2 below), we establish several 
pointwise, uniform estimates which show that the first- and the second-order 
approximations approach closely the shock wave solutions of (1.3) with suf- 
ficiently small strength. See Theorem 3.2 (for Burgers equation). Theorem 
4.2 (general conservation laws), and Theorem 5.1 (generalization to second- 
order approximation). In the last section of the paper, we discuss whether 
our results are expected to generalize to higher-order approximations. 



2. Formal Chapman-Enskog expansions 

2.1. Expanding only. In this section we will discuss two variants to 
derive a formal Chapman-Enskog expansion for (1.3), at any order. We 
begin by plugging the expansion v = "^kLo^'^^k into (1.3) while keeping u 



(1.4) 
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fixed. We obtain 

oo 

dtu + ^e^ d^Vk = 0, 

fc=0 

k=0 k=0 

The second identity above yields 

f{u) = vo, 

dfVo + dxu = -vi, 
dfVk = -Vk+i, k>l, 

which determines vq = f{u) and, for k > 1, = {—\)^d^~^{dtf{u) + 
a? dxu) , while the function u is found to satisfy 

oo 

(2.1) dtu + d,f{u) = -d, (-e)' {dtf{u) + d,u) . 

k=l 

For instance, to first order we find 

(2.2) dtu + dxf{u) = edx {dtf{u) + d^u) , 
and to second order 

(2.3) dtu + dxf{u) = edx {dtf{u) + d^u) - e^dxt{dtf{u) + d^u) . 

The corresponding traveling wave equation satisfied by solutions of the 
form 

u{x,t) = u{C), C:={x-\t)/e 

read 

oo 

(2.4) -Xu' + f{u)' = 5^ A'^-i ((-A f{u) + a") u') . 

k=l 

To first order the traveling wave equation is 

(2.5) -\u' + J{u)' = {{-Xf{u) + a") u')' 
and to second order 

(2.6) -\u' + f{u)' = ((-A f{u) + a^) u')' + A((-A f{u) + a^) u')" . 

2.2. Expanding both and v^. One can also expand both and v^, as 
follows: 

oo 

u^ = uq + eui + ... = no + ^ Uk, 

k=l 

oo 

Ve = vo + evi + ... = vq + 'Y v^. 

k=l 
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The solution at -order is defined by 

(2.7) Uk ■■= uo + eui + ... + e'' Uk- 
We also set 

(2.8) Vk ■■= vo + €vi + ... + e'' Vk. 
To first order, one can write (1.3) as 

dtuo + e dtui + d.^vo + e d^vi + O(e^) = 0, 

dtVQ + edtvi + {da^uo + ed^ui) +0(e^) = -(/(tto) + ef{uo)ui -VQ-evi) + 0(e), 
which yields the following equations: 

f{uo) -vo = 

dfUo + d^vo = 0, 

dtUi + dxVi = 0, 

dfVo + 9a; "0 = f'iuo) Ui-Vi. 

Thus 

9t«o + dxf{uo) = 0, 

Sfui + dx{f{uo) ui - dtVQ - d^u^) = 0. 
Therefore, the first-order, Chapman-Enskog expansion leads us to 
dtUi + dxf{u\) = dt{uo + eui) + dx{f{uQ) + ef'{uo)ui) 

(2.9) = e {dxtVQ + a^dxxUo) 

= e [a^dxxUQ - dttuo). 

Using that dfUo = —dxf{uo) we get 

dtui + da,f{ui) = ed4{-f'{uof + a'')dxuo)+0{e^). 

Neglecting the terms in O(e^) we may consider that ui = uq + eui is & 
solution of 

(2.10) dtui + dxf{ui) = edx{{-f'iuif + a^)dxui). 

By a similar, but more tedious calculation we can also derive the diffusive 
equation at second-order. Using (2.7) and (2.9), we have 

dtU2 + dxf{u2) = dtUi + e^dtU2 + dx{f{ui) + f{ui) U2) 

= ^{a^d'xxUQ - dttUo) + e^{dtU2 + d^^f {ui) U2)) . 

But, the second order expansion in (1.3) gives 

dtU2 + dxV2 = 0, 

dfVi + c? d^ui = /'(Si) U2-V2, 

and we get 

dtU2 + dxf{u2) = e{a^dxxUo - duuo) + dx{f {ui) U2 - V2)) 
= e[a^dxxUQ - dttUo) + e^dx{dtvi + 0? d^ui) 
= e[a^dxxUQ - dttUo) + ( - duui + dxxUi) 
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Finally, since ui = uq + eui we conclude that, to second order, 

(2.11) dtU2 + d^f{u2) = e{a^d^^ui - duUi) . 

In exactly the same manner we have, for n > 1, 
so that 

(2.12) dtUn + d^fiun) = e {a^d^^un - dttun) + 0(e"+^). 

In general, the n^'^-order equation is obtained by replacing duUn-i by deriva- 
tives with respect to x to obtain an equation of the form 

n 

(2.13) dtUn + dJ{un) = Y,^^ 

) dx ) • ■ • ) dx ^1 ) ■ 

k=l 

We will refer to this expansion as the Chapman- Enskog expansion to n^^ 
order. 

So let us for instance derive in this fashion the second order equation 
satisfied by U2- We have first 

(9 izL*?**^! = dt{dtui) = dt {-dxf(ui) + edxiia^ - f'{uif)dxUi)) 
^ ^ = -dx{nm)dtui) + edxtiia^ - nm))dxui) 

Then setting 

g[{u) =a^- f{uf and = " = g[iu)f'iu), 

we get 

dttui = -dx {f\ui){-f'{ui)dxUi + edxxgi{ui))) + edxxtgi{ui) 

= dxif {uifdxUi) - edx{f{ui) OxxQiiui)) + edxx{g'i{ui)dtui) + O(e^) 

= dx{f'{uifdxUi) - edx{f'{ui)dxxgi{ui)) + edxx{g'i{ui){-f'{ui)dxUi) +0{e^) 

= dx{f'{uifdxUi) - edxW {ui)dxxgi{ui)) - edxxxg2{ui) + 0{e^). 

Finally, since ui = U2 + O(e^), from (2.11) we obtain 

(29JS^ + dxf{u2) = edxxgi{u2) + Qx {!' ^2)dxxg\^2)) + 9^x52(^2)) • 

Setting u = U2, we can rewrite the last equation in the form 

+ ^{u)x =e((a^ - f'(uf) Ux)x 

(2.16) +e^fXu){ia^-f'{uf)ux)x), 

+ €^{0" - f'{uf)f'{u) Ux)xx. 

For later reference we record here the traveling wave equation associated 
with ([236]) 

-Xu' + f{uy={{a'-nuf)u')' 

(2.17) +^f'^u){(a'-nuf)u')')' 

+ {{a^-f{uf)f{u)u')". 
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We arrive at the main issue in this paper : Does the solution Un of (2.13) 
converge to some Umit u when n ^ oo and, if so, does this hmit satisfy the 
equation 

dtu + dxfiu) = e{a^dxxu - duu) . 
In other word, is this limit u a solution of the relaxation model (1.3) ? To 
make such a claim rigorous one would need to specify in which topology the 
limit is taken. As we are interested in the regime where shocks are present 
the convergence in the sense of distributions should be used. We will not 
address this problem at this level of general solutions, but will investigate 
the important situation of traveling wave solutions, at least as far as first- 
and second-order approximations are concerned. 

3. Burgers equation : Validity of the first-order equations 

We begin, in this section, with the simplest flux function f{u) = u'^/2. 
Modulo some rescaling x — > x — Xt/e, the traveling wave solutions u = u{x), 
V = v{x) of (1.3) are given by 

-A u' + v' = 0, 



(3.1) 



-A v' + u' 



2 



V, 



where A represents the wave speed. Searching for solutions connecting left- 
hand states M_ and f_ := f{u^) to right-hand states u+ and v+ := f{u+) 
(so both at equilibrium), we see that 

A (i(+ — U-) = v+ — V-, 

so that the component u is a solution of the single first-order equation 

(a^ — A^) u' = -{u — u-)(u — u+). 

The shock speed is also given by A = (-u^ -|- ■U-)/2. Finally, an easy calcu- 
lation based on (3.1) yields the following explicit formula for the solution, 
say u = u^{x) of (3.1) connecting «_ to It exists if and only if ?x_ > «_|_ 
and then 

(3.2) u,{x):=u. (^--^+) 



1 I / '11+ 

It will be useful to introduce the following one-parameter family of func- 
tions 

(3.3) V ^^^\{«'>' 

1 + exp^- 

in which is a parameter, not necessarily related to the speed A. Clearly, 
we have 

Note that we have for all ^ < a^, and x G M, 

U+ < (Pfj,(x) < U- 
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The following estimate in terms of the strength d := {u- — tt+) is easily 
derived from (3.3): 

Lemma 3.1. Given a > and < h < there exist constants c, C > 
such that for all fii,H2 € + h,a'^ — h) and for all x eM. we have 

(3.4) \ip^,{x)-ipi,,{x)\<CS''\x\\fii-fi2\e-'^\^\^. 

Proof. We can write 



(3.5) 
- 2 



1 + exp 



x{u. 

-2(a 



exp 



sup- 

x,k M _|_ 



1 + 



exp 



X 5 



2 ■ 



Here, the super bound is taken for |fc| < — /t and x G 
Then observe that for y > we have 



exp 



y 

2{d^-k) 



(l + «^p(-2(^))" 



< 



^ exp (-2(^) 

(-2(a^+(a^-ft))) ' 



while for 77 < we have 



exp 



( 2(^^) 



< 



1 + exp 



'2(a^-fc) 



^ exp ^2(;if3fc)^ 

< exp ■ 



This establishes the desired estimate. 



□ 



We are now in position to study the traveling waves of the first-order 
equations obtained by either the approaches in Subsections 2.1 and 2.2: 



and 



-An'+(^y=((a^-n^)^')'' 



respectively. Note that they only differ by the diffusion coefficients in the 
right-hand sides. After integration, calling Vi and Wi the corresponding 
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traveling wave solutions, we get 

(3.6) (a2 - A Vi) V; = ^ (Vi - u^) {Vi - u+) 
and 

(3.7) (a2 - Wl) W[ = ]^{Wi-u.) {Wi - n+), 

respectively. For uniqueness, since the traveling waves are invariant by trans- 
lation, we assume in addition that for example 

(3.8) uM = vm = wm = !^^±^. 

To compare the first-order diffusive traveling waves Wi and Vi with the 
relaxation traveling wave n*, we rely on monotonicity arguments. It is 
clear that the traveling waves are monotone, with Vl , W[ < and > 
Vi{x) , Wi{x) > u+, so that setting 

r_ = min u'^ — b6, r_|_ = max + b6, 

[«_l_,n_] [ti_|_,ii_] 

where 6 > is a sufficiently small constant such that r+ < a^, we find 
(a^ - r_) Wi<l {Wi - n_) {Wi - u+), 

(3.9) 

(a2 - r+) W[>^ {Wi - n_) {Wi - u+). 
Therefore, setting 

u = Wi- ipr_ , 

after some calculation we find 

1 - exp (—■nT^rzf^ 
(3.10) 2{a^ -r.)u' -u^ + u6 ) ^" <0 



1 -|- exp 



xS 

'2(a^-r_) 



We have m(±oo) = 0. As x — > ±oo the last coefficient in ()3.10p approaches 
±1 and the function u satisfies 

cu' ±u6 + B..O.T. < 0. 

So, u decreases exponentially at infinity while keeping a constant sign, and 
we deduce that u{x) 7^ for \x\ > M, for some sufficiently large M. 

Now, if u vanishes at some point xq then, thanks to the inequality (jS.lOp . 
we deduce that u'{xq) < 0. This implies that there is at most one point, and 
thus exactly one point where u vanishes, which is by (3.8) xq = 0. Therefore, 
we have sgn(x) u{x) < 0. 

A similar analysis applies to the function Wi — ipr^ and we obtain 

(3.11) sgn{x) ifTj^ix) < sgn{x)Wi{x) < sgn{x) ipr_{x), x € M. 
Concerning the function Vi, by defining 

A_ := min u, A+ := max u, 
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and 

A_ := min(A A_ , A A_|_) — bS, A_|_ := max(A A_ , A A_|_) + b5, 

where 6 > is a sufficiently small constant such that A_|_ < , we obtain in 
the same manner as above 

(3.12) sgn{x) ipA+ix) < sgn(x)yi(x) < sgn(x) 99a_ (x), x € M. 

Note that, for the same reasons, the function u = satisfies also (3.11) and 

(3.12) . Finally, since |r+ - r_l, |A+ - A_| < CS, we can combine (3.11) 
and (3.12) with Lemma 3.1 and conclude: 

Theorem 3.2. Given two reals a > M > 0, there are constants c, C > 

so that the following property holds for all u^,u^ G [—M, M] . The uniform 
distance between the traveling wave of the relaxation model and the ones of 
the first-order diffusive equations derived in Section 2 is of cubic order, in 
the sense that 

(3.13) \Vi{x) - n*(x)|, \Wi{x) - u^{x)\ <C5^ \x\ e^^'^l^l, x G R. 

Note that the estimate is cubic on any compact set but is solely quadratic 
in the uniform norm on the real line: 

(3.14) 11^1 -^^*||loc(k), \\Wi-u4loo(^K) < C'd". 

4. Validity of the first-order expansions 

We extend the result in Section 3 to general, strictly convex flux-functions. 
It is well-known that a traveling wave connecting u- to Uj^ must satisfy the 
condition u- > u+ which we assume from now on. 

Set 

(4.1) P{u) = f{u) - f{u-) - A (n - n_), 

and denote by n* the solution of the relaxation equation and by V\ and W\ 
the first-order traveling waves corresponding to equation (2.2) (i.e, (2.5)) and 
to (2.10) respectively. We have 

(a^-A^X = P(n,), 

(4.2) (a2-A/'(Vi))y/ = P(Vi), 

{a^-f'{W^f)W[ = P{W^), 

together with the boundary conditions 

\\m.u^{x) = limVi(x) = lim Wi(x) = u±. 

±CX) ±oo ±oo 

The existence of solutions to these first-order O.D.E.'s can easily be 
checked, for instance using the following implicit formula: 

Fk{u{x)) - Fk{u{0)) =x, X eR, k = 0,1,2, 
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where 



(4.3) 



Fiiu) :- 











-f{u- 


) — X{u — u-) ' 






A/'(^)) 




-f{u- 


) — A (tt — U-) ' 









To ensure uniqueness, we can impose, for example, 

U- + u+ 



u G 



u G 



u G 



(4.4) 



w*(0) = Vi{0) = Wi{0) 



Now, as was done for Burgers' equation, let us define auxilliary functions 
as the solutions of 

(4.5) {a^-fi)ip'^ = P{^^), 

with the same boundary conditions as above. For fj, < we immediately 
have 

u+ < <fij,{x) < U-, X G M. 
Setting S := {u- — u+) we get: 

Lemma 4.1. Suppose that f is a strictly convex flux-function and m_ > u_|_. 
Given a > and < h < there exist constants c, C > such that, for 
all /ii, /i2 £ (— + h,a^ — h) and for all x G M, 

(4.6) |v9^,(a;)-v?^,(x)| <C52|x||//i-//2|e-'=l^l^. 
Proof. Let il) be the solution of 

V^' = P(V') = /(V) - f{u-) - A(V - u.). 

We clearly have 

',^/x(a;) = V' 



Now, we can write 

\'Pni{x) - ¥5^2 (a;) I 



a? — ^ 



a" - ^1 
xil)'{k{x) x) 



— /xi — 112 



< C - fi2\\x\\P{tp{k{x) x))\ . 

Here, k{x) is some real number lying in the interval ^^^rz^, ^rz^^- 
On the other hand we have 

\Pi^{x))\ < CS\^|J{x) -u-\ < CS^. 

This implies that 

(4.7) YP\xx{x) - H>ti.2{x)\ < C\fii - H2\S'^\x\. 
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The behavior at ±00 is described by 

X — > +00 

and 

il^ix) ~ A;_e(-^'("-)-^) X -00. 

Since the coefficient k{x) is bounded away from and f'{u-^) — A = c+ (5 and 
f'{u^) — A = c_ (5 with C-|_ < and c_ > (bounded away from zero since 
/ is strictly convex), this completes the proof. □ 

Consider now the functions Vi and Wi the solutions of (4.2). Then, 
we have: 

Theorem 4.2. Let f be a strictly convex flux-function, M > and a > 
such that (1-2) holds in [— M, M]. Then there exist constants c, C > so 
that the following inequality holds for all u-,u+ G [—M,M] with u- > u+; 
for all X eR 

(4.8) \Vi{x) - u^{x)\, \Wi{x) - u4x)\ <C6^ Ixje-'^^l^l. 
The proof relies on the following lemma: 

Lemma 4.3. Suppose that f is a strictly convex flux-function andu- > u+. 

Assume that zj^ and z_ are the solutions of 

z'+ = R+{z+), z'_ = R.{z^), z+{0) = z^{0), 
where R+ = R+{u) and i?_ = R-{u) are any smooth functions satisfying 

(4.9) R+{u) < R^{u) <0 for all ue{u+,U-). 

Then, the two corresponding curve solutions cross at x = only, and 

z+ > Z- for X < 0, 

(4.10) 

z+ < Z- for X > U. 
Proof. If there is xq such that 2+(xo) = Z-{xo) then thanks to (4.9), 

z+{xo) < z'_{xo). 

This implies that there cannot be more than one intersection point. So, 
(0, 2;+(0)) is the only interaction point of the two trajectories, and (4.10) 
follows as well. □ 

Proof of Theorem 4.2. Setting 

A_ = min /'(u), A+ = max /'(u) 

IU+,U-] 

and 

A_ = min(A A_, A A+) — bd, A+ = max(A A_, A A+) -\-bS, 
where, 6 > is a sufficiently small constant such that A+ < a^, we have 
A_ < A /'(«)< A+ and A_ < A^ < A+, 
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and thus 

(4.11) < - A+ < - A f'{u), - A^ < - A_. 
Applying Lemma 4.3 we deduce that 

(PA. < U^, Vi, < (pA+ X <0, 
< U*, Vi, < ipA_ X > 0. 

Now, concerning the third equation in (4.2), wc set 

r_ = min f'{uf' — h5 and r_|_ = max f'{u)^ + b5, 

["+>"-] 

where 6 > is sufficiently small such that < a?. We obtain 

(4.12) < - r+ < - f'{uf, - \^ < - r_ 
and, by Lemma 4.3, 

ipT_ < u*, Wi, < (^r+ X <0, 
'Pr+ < u*i Wi, < ipr_ X > 0. 

Finally, since |A_(_ — A_|, |r-|_ — r_| < C 5, by applying Lemma 4.1, we 
obtain (4.8). This completes the proof of Theorem 4.2. □ 

5. Validity of a second-order expansion 

Our next objective is to extend the estimate in Theorem 4.2 to the second- 
order equation obtained in Subsection 2.1. 

We consider the equation (2.6) after integrating it once. The traveling 
wave connects n_ to with n_ > it+, and is given by 

(5.1) P{u) := (-A f'{u) + a2) u' + X ((-A f'{u) + a^) u')' . 
Defining first- and second- order ODE operators: 

Qiu = (a^ - Xf'{u))u' 

and 

Q2U = (a^ - Xf'{u))u' + A((a2 - A/'(ti))ti')' = Qm + X{Qiu)' . 

The solution n = V2 of (2.6) under consideration satisfies 

(5.2) Q2V2 = P{V2) 

Theorem 5.1. Let / : R — > R 6e a strictly convex flux-function and M > 0. 
Then there exist constants C, c, cq > so that the following property holds. 

For any U-,u^ € [—M, M] with n_ > U-^- and < 5 = n_ — u+ < cq, there 
exists a traveling wave V2 = V2{y) of (5.2) connecting U- to Moreover, 
this traveling wave approaches the relaxation traveling wave "u* to fourth- 
order in the shock strength, precisely: 

(5.3) \V2{x)-u^{x)\<C5^\x\e-''\''\\ x G R. 
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The estimate is only cubic in the uniform norm on the whole real line: 
(5.4) \\V2-u4loo(^)<C'5\ 
Proof. Setting 

d, = ^, and lx = d,. = ^^, 

then = ipx2 satisfies 

Qiu, = + 7a(A - f'iu,))) = PK)(1 - 7xP'iu*)), 

and a simple calculation gives 

Q2U, = PK)(l-72(/"(n,)PK)+(/'K)-A)2)) = P(n,)(l-7f(PP')'K))- 

In the same manner, the function (/j^, that is the solution of (4.5) satisfies 
the following equation 

Q.if^ = P{^^) (1 + + d^{X - fiip^))) = P{ip^) (1 + - d^P'iip^)) , 

where 

_ - 
— ji 

and 

Q2^,. = P(¥';.)(i+c^(1+c/^(/'(<^m)-A))-4(/"(<^m)^(¥'m) + (/'(^m)-A)')) 

or, equivalently, 

Q^^^ = P{^^) (1 + c^(l + df,P'{^^)) - 4(PP')'(<^^))) . 

Now, since \f{v^)-X\ < Co (5 and |/"(^^)P(^^) + (/'(^^) - A)2| < CoS\ 
then for sufficiently small S there exists a positive constant C such that the 
following property holds: by choosing /j,^ and /x_ in the form 

H+ = X'^{1 + CS'^), = A^(l - C(52), 

we obtain 

Q2^i^+ = PiVn+){^ + K+{ip^^)), where K+{ip^^) > 

and 

Q2V^^_ = P(<^M-)(1 + K-{ip^,_)), where K^{ip^_) < 0. 

Consider the corresponding functions ip^^ and and let us use phase 
plane argument. The corresponding curves 

C- : </J^_ {(Pi^_,w^_ = Q\<p^_) 

satisfy 

(5.6) A/(<^^Jti;^+^ + = P(<^^J(1 + i^+(<^^J) 
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and 

(5.7) Ai(^^_)w;^_^ + w;^_ =P(¥P^_)(l + i^_(^^_)), 
where 

We claim that the curve C+ is "below" the curve 

This is true locally near the points {u-,{)) and (u+,0), as it clear by 
comparing the tangents to the curves at these points (using (4.5)). Note 
that if A = we have u = u*. We then distinguish between two cases: 

Case 1: If A > 0, suppose that the two curves issuing from (tt_, 0), meet 
for the "first" time at some point {uq^Wq) with u+ < uq < u-. Then, 
combining (5.6) and (5.7) at this point we get 

\l(uo)w^ (^(mo) - ^N)) = P{uo){K+{uo) - K_{uo)). 
This leads to a contradiction, since 

wo < 0, ^^(^o) < ^^(^o) and P{uo){K+{uo) - K^{uo)) < 0. 

Consider now the equation (5.2) and let us study in the phase plane the 
trajectory issuing from (m_, 0) at — oo. Comparing the eigenvalues we obtain 
that the tangent at this point lies between those of the reference curves C+ 
and C-. 

In the same manner as before, we obtain that this curve cannot meet C+, 
nor C_, and necessarily converges to («+,0) as y ^ +oo. 

Case 2: If A < 0, we follow the same analysis by considering the trajec- 
tory of (2.5) arriving at (fx+jO) and the "last" intersection point. 

In both cases, we obtain the existence (and uniqueness) of the solution 
of (5.2), denoted by u = V2, and also that its trajectory called C is between 
C+ and C_. 

Note that since our equations are autonomous, by choosing ti(0) = (/?^_^ (0) = 
"^//-(O) = ('i*- + ■*^+)/2, we have 

(5.8) < u < (pfj__, x>0 
and 

(5.9) <u< ip^^ , X <0. 

Indeed, from the phase plane analysis, if for some xq G M, u{xo) = ^Pfx+{x{i) 
then necessarily w{xq) > w^_^^{xo) and then u'{xo) > </?Ji_^ (a^o)- This means 
that the curves x 1— u{x) = V2{x) and x ipfj__^_{x) have only one intersec- 
tion point, that is (0, m(0)), that satisfies in addition u'{0) > ^'^_^_{0). We 
obtain in same manner that the two curves x 1— > u{x) and x 1— > (pi^_{x) 
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have only one intersection point, that is (0,tt(0)), that satisfies in addition 

Now, using the inequahties (5.8) and (5.9) that are also satisfied by u* = 
(px^ ( since /i- < < fi-\-) , we can write 

\u:,{x) - u{x)\ < |<^At+(x) - ^Pf,_{x)\ 

< |;u+ -/i_|(52|x|e-'=l^l^ 

which completes the proof of Theorem 5.1. □ 

6. Conclusions 

For the general expansion derived in Subsection 2.2 we now establish an 
identity which connects the relaxation equation with its Chapman-Enskog 
expansion at any order of accuracy. By defining the ODE operator 



(6.1) Qnu := J2 ((-A/'(n) + a^WY 



,2^„./^(fc-l) 
k=i 

we have: 

Theorem 6.1. The traveling wave of the relaxation model satisfies 

QnU^ = Piu^:) (1 - 7^ Rniu^)) , 
where ^\ := A/(a^ — A^), and the remainders Rn are defined by induction : 

Ri := P', Rn+i := (P Rn)' for n > 1. 
Proof. Note that the ODE operators Qn satisfy 

Qn+lU = QlU + XiQnuY- 

Now, assume that 

QnU^ = P{u^) (1 - 7a Rniu^)) , 

then 

Qn+iu, = P{u,) (1 - 7A P'{u,))+X {P' (1 - 75J Rn{u.)) - P{u,)^l K{u*)) <. 
But since u'^ = ^r^^ it follows that 

Qn+lU, = P{U,) (1 -7^' {PRnYiu*)) = P{U,) {l - Rn+l{u,)) , 

which completes the proof. □ 

Theorem 6.1 provides some indication that, by taking into account more 
and more terms in the Chapman-Enskog expansion, the approximating trav- 
eling wave should approach the traveling wave equation of the relaxation 
equation (1.3). For n large but fixed it is conceivable that, denoting Vn the 
solution of QnU = P{u), 

(6.2) ||F„-«*||ioc(M) <C„(5"+^ 
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However, one may not be able to let n ^ oo while keeping 6 fixed. In fact, 
numerical experiments (with Burgers flux) have revealed that the remainders 
satisfy only 

where the constants grow exponentially and cannot be compensated by 
the factor 7". One can also easily check, directly from the definitions, that 

||i4K)||Loc <CS^n\. 

In conclusion, although we successfully established uniform error esti- 
mates for first- and second-order models, it is an open problem whether 
such estimates should still be valid for higher-order approximations. Theo- 
rem 6.1 indicates that the convergence might hold but, probably, in a weaker 
topology. 
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